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Abstract 

We investigate the 0{as) correction to e+e^ — )• J/ip + gc2 in the NRQCD factorization 
approach. A detailed comparative study between e"'"e~ J/iIj + r/c2 and e^e~ J /iIj + x'ci 
B factory energy is also carried out. After incorporating the 0{as) correction, we predict the 
cross section for the former process to be around 0.3 fb, while that of the latter about 6 times 
greater. The outgoing J/ip is found to be dominantly transversely-polarized in the former process, 
while longitudinally-polarized in the latter. These features may provide valuable guidance for 
the future experiment to examine the ^Pi or charmonium option of the X(3872) meson 
through the exclusive double-charmonium production processes. The observation potential of 
e^e~ —7- J/ip + x'ci looks bright for the current data sample of the Belle experiment, provided 
that the x'ci is indeed the narrow X(3872) state. In the appendix, we also identify the coefficients 
of the double logarithms of form ln^(s/m^) associated with all the relevant next-to-leading 
order Feynman diagrams, for the helicity-suppressed double-charmonium production channels 
e+e" J/V' + ^?c2(??c,XcO,l,2)• 
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I. INTRODUCTION 



Since a number of double charmonium production processes were first discovered at B 
factory a decade ago this topic has spurred a widespread interest Firstly, this 

production environment provides a unique and powerful means to search for the new C- 
even charmonium states, especially those X, Y, Z states, by fitting the recoil mass spectrum 
against the J/ip [ip')- The most famous examples are the discovery of the X(3940) [5] and 
the X(4160) ^] by this way. 

Another appealing reason to study double charmonium production is that it provides 
a new stage to sharpen our understanding toward perturbative QCD, especially toward 
the application of the light-cone approach 0, [sl and the nonrelativistic QCD (NRQCD) 
factorization approach (9| to hard exclusive reactions involving heavy quarkonium. 

The most famous double-charmonium production process is perhaps e"'"e^ —> J/ip + rjc- 
The original lowest order (LO) NRQCD predictions to this process [lO, lH is about one order 
of magnitude smaller than the measurement This acute discrepancy has triggered a great 



amount of theoretical investigations in both NRQCD and light-cone approaches |12|-|22|. One 



crucial element in alleviating the discrepancy between the NRQCD prediction and the data 
is the significant and positive next-to-leading order (NLO) perturbative corrections 16|, 111 . 
By contrast, owing to some long-standing theoretical obstacles, the NLO correction to this 
helicity-suppressed process in the light-cone approach has never been successfully worked 
out. As a consequence, despite some shortcomings, the NRQCD approach seems to be the 
only viable method which is based on the first principle and also systematically improvable. 

Recently, the 0{as) corrections to the double-charmonium production processes e"'"e~ — >• 
J /ijj + Xco,i,2 have also been investigated in NRQCD approach |23|, [2^]. For some production 
channels, the effect of the NLO perturbative corrections can be important. 

In this work, we plan to carry out a comprehensive study of the NLO perturbative 
corrections to the processes e'^e" — )■ J/ip + ric2 and e~^e~ — )■ J/i' + x'ciy again in the NRQCD 
factorization framework. The calculation for the former process is new, while that for the 
latter can be readily adapted from Ref. [l^]. This work should be considered as a sequel of 
Refs. HQ- 



Our interest in conducting such a comparative study is strongly motivated by the con- 
troversy recently raised at the quantum number of the X(3872) meson, whether it being 
l^^ or 2 ^ [25[. Consequently, the canonical charmonium options for the X(3872) would be 
the x'ci and r]c2, respectively. In literature, there have already existed a number of studies 
about critically distinguishing the properties of x'ci and rjc2, such as the radiative transitions 



the ?7c2(Xci) hadroproduction rates [29|], or ric2{x'ci) produc- 



tion rates in B decay |30|. By confronting the established properties of X(3872), all these 
studies tend to disfavor the ^D2 assignment of the X(3872) meson. We hope that double- 
charmonium production can also be added to the above list as a valuable means to help 
clarify the situation. 

It turns out that the production rate of e~^e~ J/ip + x'ci is about 6-7 times greater 
than that of e~^e~ — )■ J /ip -\- ric2- Based on the 1 ab~^ data currently accumulated at the 
Belle experiment, it looks promising to observe the former process, if the x'ci can indeed 
be regarded as the X(3872) particle with very narrow width. This seems to constitute the 
strong enough incentive for experimentalists to make an updated analysis of the double 
charmonium production at B factory. 

The rest of the paper is organized as follows. In Section |Tll we specify the helicity selec- 
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tion rule suited for the hard exclusive reaction e"'"e~ J /ip + ric2, and give the definition 
for the dimensionless, reduced helicity amplitudes. In Section llllt we present the leading 
order expressions for all the independent helicity amplitudes in the NRQCD factorization 
framework. In Section [TV] we first review some key technical issues about the 0{as) calcula- 
tion, then present the asymptotic expressions for the NLO perturbative corrections to all the 
encountered helicity amplitudes. The pattern of the double-logarithmic scaling violation is 
confirmed once again. In Section |Vl a comparative study is performed for both unpolarized 
and polarized cross sections between the processes e~e~^ — )■ J/ip + r]c2 and e~e^ — )■ J/ip + x'ci- 
This study may shed some light on unveiling the quantum number of the X(3872) meson in 
the future double-charmonium production experiments. Finally we summarize in Section IVTl 
In Appendix \X\ we tabulate the coefficients of the double logarithm ln^(s/m^) associated 
with all the relevant NLO Feynman diagrams, for all the double-charmonium production 
channels we have studied so far, i.e., e~^e~ J/ip + ric2{,'Uci Xco,i,2)- 



II. HELICITY SELECTION RULE AND REDUCED HELICITY AMPLITUDES 

It is often desirable to glean more information than simply present the unpolarized cross 
section for a hard exclusive reaction, especially for the double-charmonium production pro- 
cess considered in this work. It is of particular advantage to making explicit predictions 
to various J /il^ + ?7c2 production rates for different helicity configurations. The underlying 
reasons of carrying out such detailed studies are two-fold. On the experimental ground, once 
the sufficient statistics is achieved, the helicity amplitudes themselves in principle can be 
measured by studying the angular distributions of these charmonia and their decay products; 
from the theoretical perspective, it is also instructive to stay with the helicity amplitudes. 
The reason is that, for a hard exclusive reaction, the relative importance of the polarized 
cross section in a given helicity channel is dictated by the celebrated helicity selection rule 

(HSR) tail]. 

We are interested in the hard-scattering limit a/s ^ rric^ Aqcd, where a/s stands for 
the center-of-mass energy of the e^e~ collider, ttlc for the charm quark mass, and Aqcd for 
the intrinsic QCD scale. In this limit, the asymptotic behavior of the production rate for 
J /ip + r]c2 in a definite helicity configuration follows from the HSR (lol |: 



a[e+e- J/^/'lAi) + r/.alAa)] 



2N 2+IA1+A2I 



- > (1) 



cr[e+e — 7- ] 

where Ai, A2 represent the helicities carried by the J/ip, ?7c2, respectively, v denotes the 
characteristic velocity of charm quark inside a charmonium. Equation ([T]) implies that the 
helicity state which exhibits the slowest asymptotic decrease, thus constitutes the "leading- 
twist" contribution, i.e., a ~ is (Ai, A2) = (0,0). 

We have chosen by default to work in the e^e~ center-of-mass frame. Let |P| signify the 
magnitude of the momentum carried by the J/ip {r]c2), and 6 denote the angle between the 
moving directions of the J/ip and the e~ beam. Following the steps elaborated in the 
differential rate for polarized J /ip + ?7c2 production in e^e~ annihilation can be expressed as 

da[e+e- J /ipiXi) + r]c2{\2)] 



dcos6 



a 
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where ^Ai,A2 is the hehcity amphtude associated with the virtual photon decay into J /tp+rjc^ 
carrying the hehcity component (Ai, A2). 

Parity invariance can be invoked to reduce the number of independent hehcity amphtudes: 

^Ai,A2 = -^-Ai -A2, (3) 

hence the two hehcity amphtudes related by flipping the helicities of two chamonia bear the 
equal magnitude. An immediate consequence of ([3]) is that, the virtual photon decay into 
the longitudinally-polarized J/ip and ric2 is strictly forbidden by parity invariance. 

Integrating ([2]) over the polar angle 9 and including all the allowed hehcity states, it is 
then straightforward to obtain the unpolarized cross section: 



c,[e+e- ^ J/ij + = ^^^^ (2 lAil' + 2 |^i,o|' + 2 + 2 1^,2!') , (4) 

In conformity with the constraint |Ai — A2I < 1, as demanded by angular momentum con- 
servation, there are totally 4 independent hehcity amplitudes for 7* — )■ J/ip + r/c2 (Recall 
that ^0.0 = owing to parity invariance). In equation (jlj), we have also retained a factor of 
2 explicitly to account for the contributions from those helicity-fiipped states. We have also 
adopted the approximation ^ ^1 — 4r by assuming Mjj^ ^ M^^j ^ 2m c- 

In the NRQCD factorization framework, the product of two nonperturbative factors, i.e., 
the (second derivative of) wave functions at the origin for the charmonia J/ if), and ric2'. 
Rj/ipi^): and R'^^^{0), ubiquitously enters every hehcity amplitude, thereby it appears con- 
venient to define a reduced dimensionless hehcity amplitude, of which these nonperturbative 
factors are pulled out. We introduce the reduced hehcity amplitude, aAi,A2; which is related 
to the standard hehcity amplitude ^Ai,A2 as follows: 

^ 2«^/5e.eaA,,(0)fi;J0) ^ 

Ssmr. 



where CcC = |e is the electric charge of the charm quark, r = 4m^/s signifies a dimensionless 
mass ratio. To make the HSR manifest, we have explicitly stripped off a factor r^^l^^+^^^l"^) 
in ([5]), so that the dimensionless hehcity amplitude aAi,A2 is expected to scale as C(r°). 

Plugging ([5]) back into (jlj), we obtain the NRQCD predictions to the polarized production 
rate for the .J /tp{\i) + ?7c2(A2) state: 

f7[e+e-^ J/^(Ai)+r7,2(A2)] 

^ '-^^M^ ^^/.(0X2(0) r'^^^^^'-^ (1 - 4r)^ |a.,./ . (6) 

III. LO PREDICTIONS FOR THE HELICITY AMPLITUDES 

The reduced hehcity amplitude aAi,A2 can be viewed as the NRQCD short-distance coeffi- 
cient, which encodes the contribution solely stemming from the momentum region between 
iric and a/s, so can be computed reliably in perturbation theory. It is convenient to param- 
eterize it as 

«Ai,A2 - %i,A2 + — «Ai,A2- I'J 
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Our central task in this work is to decipher the 0{as) correction to the reduced amphtude, 



Ai,A2 

First we recapitulate the LO calculation. To proceed, it it most convenient to consider the 



quark amplitude 7* — > cc(^S'| ) + 00(^/^2 ) using the covariant projection technique jlOl. l32 
At LO in as, there are only four Feynman diagrams that contribute, one of which is depicted 
in Fig. [lb) (For simplicity, we have neglected the QED fragmentation diagrams, whose effect 
appears to be modest). 

It is straightforward to project out 4 corresponding helicity amplitudes from the decay 
amplitude for 7* — )■ J/ip + ric2- One then follows equation ([5]) to read off each of the LO 
reduced helicity amplitudes: 



(0) _ r±(l-4r)f (Ai,A2) = (±l,0) 



\ (Ai, A2) = others. 



(8) 



For some accidental reason, only the J/iIj{±1) +77c2(0) channel has the non- vanishing cross 
section at LO in as- Substituting (|8]) into (|6]), we find agreement with the QCD part of the 



LO prediction for the unpolarized cross section first given in Ref. [10 . 



IV. NLO PERTURBATIVE CORRECTIONS TO THE HELICITY AMPLITUDES 



We start this section by first sketching some technical issues about the NLO perturbative 
calculations, followed by presenting the asymptotic expressions of the 0{as) corrections to 
all the reduced helicity amplitudes. 
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A. Outline of the calculation 



At NLO in as, there are 20 two-point, 20 three-point, 18 four-point, and 6 five-point one- 
loop diagrams for the process 7* — )■ cc{^S[^^) + cc{^D^^), some of which have been illustrated 
in Fig. [U The calculation is quite similar to our preceding works on double charmonium 



exclusive production, i.e., e+e~ J/ip + Xco,i,2 124] and the e~^e~ — )■ J/tp + rjc |2l|, so here 
we will only present a very brief description. 

We adopt dimensional regularization to regularize both UV and IR singularities. We 



follow the 't Hooft-Veltman prescription for 75, as detailed in Ref. [33|. In projecting out the 



D-wave orbital angular momentum state, one needs expand the relative quark momentum 
q to the quadratic order. Our strategy is to follow the method of region js^ via directly 
deducing the NRQCD short- distance coefficients, rather than resorting to the much more 
expensive matching calculation, through making the expansion in q before carrying out the 
loop integration. 

The only nontrivial technical problem is that one may encounter some unusual one-loop 
integrals which generally contain the propagators of cubic power, as a consequence of taking 
the derivative over q twice. The Mathematica package FIRE j35| and the code Apart ^ 
are utilized to reduce these unconventional higher-point one-loop tensor integrals into a set of 
masters integrals. With the aid of the integration-by-part algorithm and partial fractioning 
technique, it turns out that all the encountered master integrals are nothing but the standard 
2-point and 3-point one-loop scalar integrals, whose analytic expressions have already been 



tabulated in Ref. 17 



When adding the contributions of all the diagrams, and after renormalizing the charm 
quark mass and the QCD coupling constant, we finally end up with both UV and IR finite 
NLO expressions for the decay amplitude 7* — )■ cc{^S[^^) + cc{^D^^). Since everything 
becomes finite, we can safely return to the 4 spacetime dimensions, and readily project out 
all the required hehcity amphtudes. 

In Ref. js^], an all-order-in-as proof for exclusive quarkonium production has been out- 
lined in the NRQCD factorization context. It argues that at lowest order in v and to all or- 
ders in as, NRQCD factorization holds for the exclusive production of a S'-wave quarkonium 
plus any higher-orbital-angular-momentum quarkonium in e"'"e~ annihilation. Our explicit 
calculation confirms that the NRQCD short-distance coefficients affiliated with S'-wave plus 
D-wave charmonia production are indeed IR finite at NLO in as, which is compatible with 
what is asserted in 13711. 



B. Analytic expressions of NLO helicity amplitudes 

The CL^Xi X2 complex- valued, whose analytic expressions are in general quite lengthy. 
Rather than reproducing their cumbersome-looking expressions here, we are content with 
presenting their numerical values over a wide range of r, as shown in Fig. [2l 

As a matter of fact, it seems much more illuminating to know the asymptotic behaviors 
of the helicity amplitudes in the limit i/i rric. At NLO in as, one anticipates to see the 
logarithmic scaling violation to the naive power-law HSR as indicated in ([T]). Furthermore, 
conducting the asymptotic expansion in NRQCD short-distance coefficients is theoretically 
appealing, since it is equivalent to disentangling the contributions occurring at the "hard" 
scale (virtuality ~ s) from the "lower-energy" coUinear/soft sectors (virtuality ~ ""^c)? 



6 



(0,1) (1,0) (1,1) 




Vs" (GeV) 

FIG. 2: Variation of the 0{as) reduced helicity amplitudes 0,^X1 X2 ^^^^ V^- We take fi = ^ 
and rric = 1.4 GeV. The sohd curves correspond to the exact NLO results, and the dashed curves 
represent the asymptotic ones as given in The vertical mark is placed at the B factory energy 
= 10.58 GeV. 



which one can intimately link the NRQCD factorization approach and the light-cone ap- 



proach jss i^]. Such an asymptotic expansion on the NRQCD hard coefficients has been 
carried out for a number of exclusive double-quarkonium production processes, either in 
e~^e~ annihilation 21, 24, 42|, or from bottomonium decay [42'-'45'l, and some general pat- 
tern about logarithmic scaling violation has been recognized [24, 42.]. 

According to the definition in ([6]), we find the asymptotic expressions of four reduced 
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NLO helicity amplitudes for Y ^ J/4' + Vc2 to be 



"0,1 



= _ ^ j 7 ln2 r - 2(22 - 43 In 2) In r - — - Un'' + 1^^^ In 2 - 41 In^ 2 

asym 96 ^ '3 3 



+2i7r(71nr- 22 + 431n2) }, (9a) 



"1,0 



= — J 191nV + 19(l + 21n2)lnr + 12/3of In ^ + - — + — ln2 

asym 48 V 3 / 3 3 



1291n2 2 + i7r 381nr + 12/3o + 19 + 381n2 }, (9b) 



= -^|21n^r + (19-231n2)lnr-l^ + ^-ll^ln2-liln^2 
asym 24 ^ '323 2 



+OT(41nr + 19-231n2) j, (9c) 
= _^| lnV-7(2-31n2)lnr---^ + 511n2--ln22 



"1,2 



asym 24 | ' '26 2 

+Z7r(21nr- 14 + 211n2) j, (9d) 

where fi is the renormahzation scale, and (3o = ~ one-loop coefficient of the 

QCD P function, and nj = 4 denotes the number of active quark flavors. In contrast with 
(IH]), all the four reduced helicity amplitudes receive non- vanishing 0{as) corrections. 

Note that the /3oln(4/i^/s) term only resides in the (±1, 0) channel, since this is the only 
channel that has a non- vanishing tree-level amplitude. For the sake of comparison, all these 
asymptotic results of the reduced helicity amplitudes are also shown in Fig. [21 in juxtapose 
with the corresponding exact NLO results. These asymptotic results appear to converge 
with the exact ones decently well even at relatively lower ^/s, say, at B factory energy. 

From Eq. (|9]), one sees that the leading scaling violation is due to the double logarithm 
In^ r. This constitutes another example to further corroborate the earlier conjecture: the 
occurrence of In^ r in the one-loop NRQCD short- distance coefficients is always affiliated with 



the helicity-suppressed channels in exclusive double charmonium production processes [24, 



42 . 



V. PHENOMENOLOGY 

Aside from a comprehensive analysis on the process e'^e'^ J/ip + ric2 at the B factory, 
in this section we also target at a detailed investigation on the process e"'"e" J /ip + x'cii 
where the essential elements have already been set up in j23| . This study is largely motivated 
by the recent concern about the nature of the X(3872) meson, whether its quantum number 
is 1"'"+ or 2 '". The canonical charmonium options for the X(3872) are the Xci and ?7c2, 
respectively. With this consideration in mind, our study may provide some useful guidance 
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e e ->J/i^+qc2 61 




3456789 10 3456789 10 

H (GeV) 11 (GeV) 

FIG. 3: The renormalization scale dependence of the LO and NLO cross sections for e~^e~ — >• 
J/ip + r]c2 (left panel) and e~^e~ J/ip + x'ci (right panel) at -y/i = 10.58 GeV. The band is 
obtained by varying the charm quark mass in the range rric = 1.4 it 0.2 GeV. 



on unveiling the quantum number of the X(3872) through the dedicated double- charmonium 
experiment at future B factory. 

Before making concrete predictions to the production rates for e+e~ J/ip + Vc2, we 
first address a subtlety about utilizing IQ. For the J/ip{±l) + ?7c2(0) helicity state, the 



±i,o'^±i,o 



canonic way of identifying the 0{as) correction is by replacing |a-|-i^o| with 2^3f?[a 
Because of the null tree-level amplitudes for the remaining helicity configurations, the 0{as) 
correction to the total cross section is solely from the (±1, 0) channels. On the other hand, 
the polarized cross section for each helicity state is also a physical observable by itself. 
Since the non- vanishing amplitudes are first generated at 0{as) for the remaining helicity 



states, it is thereby also sensible to interpret |aAi,A2 I' = (^) Wx!,xS fo^^ (^l'^2) 7^ (±1,0). 
Though formally of order a^, these new pieces do constitute the leading contributions to 
the respective polarized cross sections, thereby still being consistent ^. We will follow this 
viewpoint in presenting our NLO predictions. 

In the numerical analysis, we take ^/s = 10.58 GeV, and the charm quark pole mass 
TTLc = 1.4 ± 0.2 GeV. The fine structure constant is chosen as a{^/s) = 1/130.9 [l9[. The 
running QCD strong coupling constant is evaluated by using the two-loop formula with 
= 0.338 GeV 16, T3|. The nonpertubative input parameters, i.e., (the derivatives of 



2 



MS 



the wave function at the origin for J/ip, XcIj Vr/i-, also suffer from a large amount of 
uncertainty. Their values have been compiled in Ref. [46|, which were deduced from several 
different potential models. We choose to use those given by the Buchmiiller-Tye potential 



model 47^: \Rj/JO)? = 0-81 GeV^ \R', (0)P = 0.102 GeV^ and \R''M? = 0-015 GeV 



Another important source of uncertainty for the NLO predictions stems from the scale 
affiliated with the strong coupling constant. As is well known, the scale ambiguity is a typical 
nuisance of NRQCD factorization approach, reflecting the fact that two disparate hard scales. 



^ The uncalculated oidei-v^ corrections to e~^e^ — > J/ip+Vc2 generally also lead to non- vanishing amplitudes 
for those helicity states other than (±1,0), which are potentially as important as the 0{as) corrections. 
For the leading contributions to the corresponding polarized cross sections, one needs also include the 
interference effect between the radiative and relativistic corrections. 
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FIG. 4: Variation of the LO and NLO cross sections for e'^e — )• J ftp + r]c2 (left panel) and 
e'^e'^ J/ip + x'ci ( right panel) with ^/s. nif. is fixed as 1.4 GeV and the renormahzation scale is 
taken as /u = 



TABLE I: Polarized and total cross sections (in unit of fb) for e^e — >• J /if) + 7/c2- We have taken 
rric =1.4 GeV. In the rightmost column, we also give the K factor for the unpolarized cross sections. 
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NLO 


0.14 
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1.6 X 10^5 


0.29 








LO 
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0.15 


1.57 






NLO 


0.12 


6.8 X 10-^ 


1.1 X 10"^ 


7.5 X lO"*^ 


0.23 





i/s and m^, are entangled together in NRQCD short-distance coefficients. In fact, the lesson 
gained from the light-cone approach strongly suggests that it is rather implausible to set the 
scales entering all as in NLO short- distance coefficient to be unanimously around y/s j42| . 
Although we are unable to circumvent the scale ambiguity problem in the confine of NRQCD 
approach, we may allow the /i to float in the range between 2mc and ^/i, hoping that the 
most trustworthy prediction may interpolate in between. 

In Fig. El we show the total cross sections both for e^e~ — )■ J /'ip+rjc2 and e'^e~ — )■ J/t/'+Xci 
as a function of /x, which also take into account the error due to the uncertainty of irtc. One 
sees clearly that, after incorporating the NLO perturbative correction, the scale dependence 
of the cross section has been reduced. In Fig. |H we also plot the LO and NLO total cross 
sections for e+e" J /ip + rjc2 as a function of y/s. Since the leading-twist (0, 0) channel is 
absent in both of the processes, the cross sections drop very rapidly (oc 1/s^) as y/s increases. 

In Tables [T] and [TTl we also list the predictions for the polarized cross sections from each 
individual helicity channel as well as the total cross sections for e^e~ J/ip + ric2 and 
e~^e~ — 7- J/ip + Xci) with rric flxed at 1.4 GeV but with the renormahzation scale chosen 
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TABLE II: Polarized and total cross sections (in unit of fb) for e^e — )• J/%1) + x'd- We have taken 
rric = 1.4 GeV. In the rightmost column, we also list the K factor for the unpolarized cross sections. 
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LO 
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NLO 


-0.0063 


0.73 


0.042 
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at three different points. The impact of the 0{as) corrections to the total cross sections 
is surprisingly alike for both double- charmonium production processes: incorporating the 
NLO pertubative correction decreases the LO cross sections by about 6% for /i = 2mc, or 
increases the LO result by roughly 30% for /i = A/i/2, or enhances the LO result by about 
60% for /i = a/s. Hence, the relative importance of the 0{as) correction increases as ix 
increases. Taking the medium value = A/i/2, our NLO predictions for the cross sections 
of e^e~ — )■ J /ip + rjc2 and e'^e~ J /ifj -\- x'ci ^-t \/i = 10.58 GeV reach about 0.3 fb and 1.9 
fb, respectively. The production rate of the latter process is about six times more copious 
than that of the former. 

Inspecting Tables [T] and [Tll one observes an interesting hierarchy of the polarized cross 
sections for the above processes. The overwhelming contribution for the J/ijj+r]c2 production 
comes solely from the (±1,0) state, with J/ip transversely polarized, while the dominant 
hehcity channel for producing J/ip + x'd is instead the (0, ±1) state, with J/ip longitudi- 
nally polarized. This characteristic of J/ip polarization may serve as a benchmark for the 
future double-charmonium production experiment to unveil the nature of the X(3872), pro- 
vided that the prospective Super B factory can observe sufficient number of X(3872) events 
recoiling against J/ip. 

We can ascertain the observation potential of these two exclusive double charmonium 
production processes. Thus far. Belle experiment has accumulated about 1000 fb^^ data. 
Taking a[e+e- J/ip + 7/^2] ^ 0.23 - 0.33 and a[e+e~ -> J/ip + x'd] ~ 1-53 - 2.16 fb 
from Tables HI and [ITl we estimate that roughly 230 — 330 J/ip + rjc2 events, and 1500 — 2200 
J + Xci events have been produced at Belle at around ^/s = 10.58 GeV. 

Since neither of the masses and decay patterns of ?7c2 and x'd known, it seems exper- 
imentally impossible to simultaneously reconstruct the J/ip and the ric2{x'd) signals. The 
viable experimental way is to only reconstruct the J/ip ^ event, and fit the recoil mass 
spectrum against the J /ifj to estimate the number of rjc2{x'd) peak events. This method 
does not depend on the concrete decay modes of ric2{x'd)i ^^^d is particularly suitable for 
the limited statistics of signal events like in our case. In fact, this method has already been 
used routinely by the Belle collaboration to search for the double charmonium production 
processes e^e~ J /if) + C-even charmonium. 
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In fitting the recoil mass spectrum, the net detection efficiency for e~^e — )■ J / i' + Vc2{x'ci) 
may reach around 4% (Similar for ?7c2 and x'd; with the reconstruction efficiency for J/ip ^ 
/+/- included Q). As a very crude estimate, the number of observed e^e — y J/'4'~^Vc2{x'ci) 
events are expected to reach (230 - 330) x 4% = 9 - 13, and (1500 - 2200) x 4% = 60 - 90, 
respectively. Since only the J/ip is reconstructed, the background level can be a little higher. 
With only around 9 — 13 observed J/ip+'r]c2 events, it appears difficult to observe a significant 
signal with current 1 ab~^ Belle data. If x'ci is indeed the X(3872) meson, thanks to the 
very narrow width of the X(3872), it seems possible to observe the 60-90 J/ip + x'ci signal 
events with the current Belle full data sample, which may provide a strong incentive for 
updating their earlier e'^e'^ — )■ J/ip+ charmonium analysis with only 673 fb~^ data [sl. |49|. 
This study may hopefully help one to better understand the mechanism about the X(3872) 
production. On the other hand, in case that x'ci is not X(3872) and its width is no longer 
narrow, it may be difficult to observe a significant signal at the moment. 

In any event, a much larger data samples seem to be called for to arrive at some definite 
conclusion. In the prospective Super B factory, which may reach an integrated luminosity of 
50 ab~^ by 2022, it seems feasible that the processes e^e~ J /if) + t]c2{x'ci) ^iii eventually 
be observed, and the quantum number of X(3872) also hopefully be nailed down. 

VI. SUMMARY 

In this work we have calculated the complete NLO perturbative corrections to e~^e~ — )■ 
J/ip + 'r]c2 within the NRQCD factorization framework, and found that the 0{as) correction 
can be sizable for a larger value of the renormalization scale. Our calculation indicates that 
the dominate contribution to the total cross section comes from the helicity states (±1,0). 

We have also carried out a comparative study for e+e~ — )■ J /ip+rjc2 and e'^e'^ — )■ J /ip+x'ci 
at B factory energy, with the specific motivation that these double charmonium production 
processes may provide useful means to unveil the quantum number of the X(3872) meson in 
the future B factories. It turns out that the production rate of the latter process is about 6-7 
times greater than the former, therefore it is possible to observe the J /ip + x'ci signals based 
on the current 1 ab^^ BELLE data sample, if the Xci is indeed the very narrow X(3872) 
particle. The dominantly transverse polarization of the J/ip is also a useful indicator for 
identifying the X(3872) with x'd- 

A necessary extension of our current work is to further incorporate the leading relativistic 
correction to the above double-charmonium production processes, whose effects might be as 
important as the radiative corrections. The study along this direction is of some theoretical 
interest, especially regarding that, until today the Oiy"^) calculation hardly exists for the 
production processes involving the P, D-wave quarkonium. 

It has been recently realized that the 0{as) NRQCD short-distance coefficients for the 
helicity-suppressed double charmonium-production processes are often plagued with the dou- 
ble logarithm of form ln^(s/m^). This symptom is likely intertwined with some long-standing 
failure of applying the light-cone approach to the "higher-twist" hard exclusive reactions 
beyond tree level. To expedite a better understanding and controlling of these double 
logarithms, we have scrutinized all the NLO diagrams for the various exclusive double- 
charmonium processes that have been studied so far, e.g., e+e~ — t- J/V' + '7c2('7c, Xco,i,2), and 
singled out those that contain the double logarithm, and enumerate their coefficients for 
each individual helicity states. No simple and general pattern has been recognized yet. 
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FIG. 5: The NLO diagrams that contain the double logarithm ln^(s/m^) (in Feynman gauge) for 
the process 7* — )■ J/V" + ?/c2 (??c; Xc0,i,2)- The respective charge-conjugated diagrams have been 
suppressed. 
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Appendix A: Anatomy of the NLO diagrams containing double logarithm in various 
double-charmonium production processes 

Although the NRQCD factorization approach and the light-cone approach are based on 
two completely different expansion strategies, they can be intimately linked for the exclusive 
double charmonium production processes in the limit s ^ m^. As exemplified by an anatomy 
of the 0{as) correction to the Be electromagnetic form factor j42|, the light-cone approach 
can be efficiently utilized to reproduce the asymptotic behavior of the 0{as) NRQCD short- 



distance coefficient through the idea of refactorization |39H41 . 

It is worth emphasizing one important limitation of the refactorization procedure. So far it 
can only be successfully applied to the "leading-twist" reactions involving quarkonium, that 
is, with the HSR-favored helicity configurations. A general characteristic about this class of 
of processes is that, at the NLO in a^, only the single collinear logarithm of In ^ arises in 
the NRQCD short- distance coefficient, which can be resummed to all orders in as with the 
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aid of tlie Brodsky-Lepage evolution equation. Tliis feature has been exphcitly verified in 



several examples, e.g., 7* r^c+T [38|, |4l|, l50| , 7* + ^^ |4l|, 7* 7/^^(0) + Xco,2(0) [2^ 

However, the majority of phenomenologically relevant double- char monia production pro- 
cesses are of "higher-twist" nature, exemplified by 7* — )■ J/4' + ric2{f]c, Xco,i,2)) for which some 
of or all of possible helicity channels are of HSR-suppressed type. Until today, it is still not 
clear how to consistently compute the NLO radiative correction to this type of processes in 
the light-cone framework, due to some long-standing problem such as the inevitable emer- 
gence of the endpoint singularity. By contrast, NRQCD factorization approach, which is 
based on the velocity expansion rather than the twist expansion, serves as the only viable 
tool to investigate the 0{as) corrections to this kind of helicity-suppressed hard exclusive 
reactions. 

An abnormal feature about the 0{as) NRQCD short-distance coefficients for the helicity- 
suppressed double charmonium production processes seems to be the frequent emergence of 
the double logarithm ln^(s/m^), which seems to seriously deteriorate the convergence of 
the perturbative expansion in NRQCD factorization ^. This double logarithm appears to 
result from the overlap between the "collinear" and "soft" regions in the NRQCD short- 
distance coefficient and seems deeply related to the afore-mentioned endpoint singularity 
problem. Although it is feasible to reproduce the closed form of the double logarithms for 
each concrete process, the lack of a thorough understanding of their behavior in higher order 
in as prevents one from systematically controlling them, i.e., to resum them to all orders in 
as- 

A useful starting point is to anatomize all the double logarithms of form ln^(s/m^) 
from the existing NLO calculations for the various exclusive double charmonium pro- 
cesses. For this purpose, we examine all the NLO Feynman diagrams for the processes 
e~^e~ — )■ J / ip + r]c2{i]c, Xco,i,2)! sort out those that contain the double logarithm, and enumer- 
ate their coefficients. 

A careful examination reveals that (almost) all the relevant NLO diagrams that contribute 
the double logarithm ln^(s/m^) (in Feynman gauge) for the above processes are included in 
Fig. O which was first identified in |42| . 

In Table IIIIl we tabulate the double logarithms that come from each individual diagram 
in Fig. O with the color-structure specified, for the processes e~^e^ J/ip + ?7c2('7c)- 

The reaction e'^e" — t- J/ip{\i) + ?7c(A2) only possesses one independent helicity config- 
uration (±1,0). Each diagram contains a non- vanishing double logarithm. Summing up 
all their contributions, and multiply by 2 to account for the contribution from the charge- 
conjugated diagrams, we recover the net double logarithm in the NLO short- distance coef- 
ficient: §|nj|a^gyjn = ^ In^ as given in equation (37) of Ref. (43 |. 

All the four polarized production channels e+e~ — t- J/ip{\i) + rjc2{^2) are helicity- 
suppressed. Adding up the contributions from all the diagrams in Fig. [5l and multiply by 
2, we recover the net double logarithm in the NLO reduced helicity amplitudes ^^^'1^2 Laym' 
as given in equations ([9]). Note that the double logarithms from the (±1, 0) channel has the 
similar diagram- by-diagram structure as its counterpart for e"'"e~ J /tp + r]c. Neverthe- 



^ In fact, one will encounter the double logarithm even in the single-quarkonium exclusive production 
process such as 7* ?7c + 7, once the corresponding 0{as) NRQCD short-distance coefficient is expanded 
to NLO in ml/ s. This can be readily checked by using the analytic 0{as) short-distance coefficient given 
in [silsBl- 
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TABLE III: The coefficients of tlie double logarithm associated with each diagram in Fig. \5\ from 
the various helicity states (Ai, A2) in the processes 7* — )• J/ip + r]c2{ric)- The two types of color 
factors are denoted by Ci = Cp and C2 = Cf{Cf — ^Ca), where Cp = | and Ca = 3. 
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h) 


c) 
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e) 
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h) 
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Color Factor 


Ci 
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C2 
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(1,0) 
(0,1) 
(1,1) 
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3 

64 


9 
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3^3 
128 

3^3 
256 

3^6 
256 


3 

256 

3^3 
256 


9 

64 


3 

128 


21 
128 

3^3 
128 


9 
128 


3 
128 


9 
128 


3 
128 


Vc 


(1,0) 


3 

64 


3 
32 


3 
128 


3 

32 


3 

64 


9 

64 


3 

64 


3 

64 


3 

64 


3 

64 




k) 



FIG. 6: The extra NLO diagram relative to Fig. [5] that also contributes the double logarithm 
ln^(s/m^) (in Feynman gauge) to the helicity channel (±1, ±2) for 7* — >• J /il^ + Xc2- The respective 
charge-conjugated diagram is not shown. 

less, for the remaining three helicity channels, a great simplification arises that many of the 
diagrams in Fig. [5] do not contain double logarithm. 

Next we examine the double logarithms that are affiliated with the polarized production 
processes e^e~ — )■ J /'ip{Xi) + Xcj{,^2) {J = 0, 1, 2). Excluding two helicity-conserved channels 
7* — )■ J/ip{0) + Xco,2(0), Table [IV] lists the structure of double logarithms in the remaining 
eight helicity-suppressed channels. A noteworthy thing is that, for the helicity channel 
(±1, ±2) of 7* — )■ J/?/' + Xc2, besides the diagrams shown in Fig. one extra NLO diagram 
as shown in Fig. [HI also need be included. The color factor associated with the diagram in 
Fig. is ^Cf, and the corresponding coefficient is — ^. 

Adding up the contributions from all the diagrams in Fig. |5] for the reactions e^e" — )■ 
J/ipi^i) + Xco,i,2(A2), together with that in Fig. |6]for the (±1, ±2) state, and multiply the 
results by 2, we then reproduce the net double logarithms in the NLO reduced helicity 
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TABLE IV: The coefficients of the double logarithm associated with each diagram in Fig. O from 
the various helicity states (Ai, A2) in the process 7* — )• J/V' + Xc0,i,2- The two types of color factors 
are represented by Ci = Cp and C2 = Cf{Cf — ^Ca)- 



Diag 


rams 


a) 




c) 


d) 


e) 


/) 


9) 


h) 





3) 


Color Factor 


Ci 


Ci 


Ci 


C2 


C2 


C2 


C2 


C2 


C2 


C2 


XcO 


(1,0) 


7 

384 


1 
192 


1 
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3 
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3 
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3 
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3 
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3 
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3 
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Xcl 
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9 
512 


3 
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9 
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3 
512 


9 
512 


9 
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3 
512 


9 
512 


3 
512 


9 
512 
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512 
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512 
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512 
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1024 
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512 


3 
512 
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512 
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3 
512 
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3 
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3 
256 


3 
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9 

128 


3 
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3 
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amplitudes, K^^ ^^{r, ^)asym given in equations (14) through (16) in Ref. [2^]. 

An obvious observation from Tables UTTl and HVl is that the specific structures of the double 
logarithms are pro cess- dependent. It seems that no simple, unified pattern can be readily 
recognized. 
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